Steinbach et at. BMC Biophysics 2013, 6:10 
http://www.bionnedcentral.conn/2046-1 682/6/1 0 



Biophysics 



RESEARCH ARTICLE Open Access 



Molecular basis of HHQ biosynthesis: molecular 
dynamics simulations, enzyme kinetic and surface 
plasmon resonance studies 

Anke Steinbach^^ Christine K Maurer^^ Elisabeth Weidel^ Claudia Henn^'^, Christian Brengel\ 
Rolf W Hartmann^'^ and Matthias Negri^' 



Abstract 

Background: PQS {Pseudomonas Quinolone Signal) and its precursor HHQ are signal nnolecules of the P. aeruginosa 
quorum sensing systenn. They explicate their role in nnannnnalian pathogenicity by binding to the receptor PqsR that 
induces virulence factor production and biofilm formation. The enzyme PqsD catalyses the biosynthesis of HHQ. 

Results: Enzyme kinetic analysis and surface plasmon resonance (SPR) biosensor experiments were used to 
determine mechanism and substrate order of the biosynthesis. Comparative analysis led to the identification of 
domains involved in functionality of PqsD. A kinetic cycle was set up and molecular dynamics (MD) simulations 
were used to study the molecular bases of the kinetics of PqsD. Trajectory analysis, pocket volume measurements, 
binding energy estimations and decompositions ensured insights into the binding mode of the substrates 
anthraniloyl-CoA and p-ketodecanoic acid. 

Conclusions: Enzyme kinetics and SPR experiments hint at a ping-pong mechanism for PqsD with ACoA as first 
substrate. Trajectory analysis of different PqsD complexes evidenced ligand-dependent induced-fit motions 
affecting the modified ACoA funnel access to the exposure of a secondary channel. A tunnel-network is formed in 
which Ser317 plays an important role by binding to both substrates. Mutagenesis experiments resulting in the 
inactive S317F mutant confirmed the importance of this residue. Two binding modes for |3-ketodecanoic acid were 
identified with distinct catalytic mechanism preferences. 



Background 

Quorum sensing (QS) is a chemical cell-to-cell commu- 
nication system in bacteria ruled by small extracellular 
signal molecules. It coordinates the social life of bacteria 
by regulating many group-related behaviours, such as 
biofilm formation and virulence factor production [1-5]. 
Anti-QS has been recognized as an attractive strategy in 
the fight against bacteria [6] based on anti-virulence and 
anti-biofilm action and not on bacterial killing. 

The opportunistic Gram-negative pathogen P, aeruginosa 
is a good model to study the complexity of QS systems 
[1,4]. At least three distinct QS pathways are known which 
regulate in a hierarchical manner the QS -dependent target 
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gene expression. The first two QS systems, las [7] and 
rhl [8], utilize A^-acyl homoserine lactones (C4- and C12- 
AHL) and the receptors LasR and RhlR [9]. The third 
QS -system is 2-alkyl-4-hydroxyquinoline (HAQ) -dependent 
and specific for P, aeruginosa and some Burkholderia 
strains [10-12]. PQS (Pseudomonas Quinolone Signal) 
and to a lesser extent its precursor HHQ (2-heptyl-4- 
hydroxyquinoline) activate PqsR [13-15]. 

A key enzyme of the PQS biosynthesis pathway is PqsD 
(PQB biosynthetic 3-oxoacyl-[acyl-carrier-protein] [ACP] 
synthase III; EC 2.3.1.180), which catalyses the formation 
of HHQ by "head-to-head" decarboxylative condensa- 
tion of anthranilate (as anthraniloyl-CoA; ACoA) and 
|3-ketodecanoic acid (|3K) [16,17]. 

Several groups have proven that P, aeruginosa pqsD 
knock-out mutant as well as PQS -deficient P, aeruginosa 
strains have an attenuated pathogenicity in nematode 
and mouse models evidencing the significance of PQS 
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signalling in mammalian pathogenesis [18]. Increased 
PQS levels have been detected in lungs of cystic fibrosis 
patients supportive for an active role of QS in chronic 
lung infections [19-21]. These findings and in particular 
the recent identification of the first class of PqsD inhibi- 
tors that reduce biofilm and virulence factor formation 
in P, aeruginosa validates PqsD as a target for the devel- 
opment of anti-infectives [22]. 

PqsD is a homodimeric bi-substrate enzyme with high 
structural similarity to FabH and other |3-ketoacyl-[ACP] 
synthases III (KAS III). They share a common thiolase 
fold (apapa), a long tunnel to the active site, and the 
same catalytic residues [23-25]. Three PDB structures of 
PqsD exist [26]: as apoform (3H76), as Cysll2-ligated 
anthranilate (CSJ) complex with ACoA molecules in the 
primary funnel (3H77) and as Cysll2Ala mutant in 
complex with anthranilic acid (3H78) [23]. In all three 
structures the catalytic centre is accessible by two chan- 
nels in L-shape: the primary CoA/ACP-funnel, and the 
shorter secondary channel (Additional file 1: Figure. SIl). 
However, the molecular details of ACoA access and, in 
particular, the binding mode and the subsequent in- 
corporation of pK are unknown. 

Knowledge of the kinetics and of the conformational 
flexibility of an enzyme can significantly contribute to a 
successful rational drug design [27-29]. Herein we study 
the molecular basis of PqsD and the HHQ biosynthesis 
combining experimental and in silico methods. Enzyme 
kinetic analysis and surface plasmon resonance (SPR) 
biosensor experiments were used to determine the 
mechanism and the substrate order of the biosynthesis; 
comparative analysis of PqsD to homologous KAS -III 
enzymes was useful to identify domains specific for 
PqsD functionality. Molecular dynamics (MD) simula- 
tions were carried out to explore the binding modes of 
ACoA and pK as well as the conformational flexibility 
of PqsD. 

Results and discussion 

Knowledge of enzyme kinetics for multi-substrate reac- 
tions is helpful to set up and interpret MD simulations. 
We performed biochemical and biophysical studies to 
determine the underlying kinetic mechanism of PqsD. 

Biochemical and biophysical characterization hint at 
ping-pong kinetic mechanism of PqsD 

Enzyme kinetic studies were performed using a 96-well 
format-based in vitro assay with the purified enzyme PqsD 
to determine the kinetic parameters for each substrate. 
Optimum enzymatic reaction conditions were determined 
in advance. Plotting product formation versus time re- 
vealed that a reaction time of 4 min in combination with 
an enzyme concentration of 0.25 \iM was suitable and that 
values in linear range within the progress curve could be 



obtained (data not shown). The initial velocity (v) was cal- 
culated by dividing the product concentration by the reac- 
tion time. Plotting the data with GraphPad Prism 5 
software resulted in an array of parallel lines in the 
Lineweaver-Burk-Plot and a common Y-axis-intercept in 
the Hanes-Woolf-Plot (Figure lA and IB). This suggests a 
ping pong kinetic mechanism for PqsD as described also 
for other Claisen condensing enzymes [30]. The results 
were in agreement with reported mass spectrometric 
[31], structural [23], and surface plasmon resonance 
(SPR) based [32] studies revealing the formation of an 
anthraniloyl-PqsD intermediate with concomitant re- 
lease of the first product CoA before binding of the sec- 
ond substrate p-ketodecanoic acid (pK). 

The /<M data (ACoA 0.875 ±0.140 (iM, pK 1300 ± 
158 \iM) correlate well with the /<d values determined 
with SPR by our group (ACoA 1.08 (iM, pK 2.95 mM) 
[32]. Also, the kinetic parameters, derived mutually 
varying both substrates (Figure IC) and fitting the data 
with the ping-pong Equation (1) are within the range of 
the apparent values determined by Pistorius et al. 
(^^M app,m< = 598.5 ± 106 (iM; V^^ax = 495.8 ± 37.5 fmol 

HHQ/min/pmol PqsD; /Ccat (PqsD as monomer) = 0.01 s'^). 

V = Vmax [^^^^ . ^ (1) [33] SPR bio- 

KmACoA+[ACoA] h + f ^ j j 

sensor assays were performed to assess the influence of 
substrate addition order on the HHQ product formation. 
Firstly, as recently reported [32], PqsD was immobilized to 
the SPR chip and ACoA injected; the increase in the 
response-line preserved also after washing was an indi- 
cator for the covalent linkage of anthraniloyl to Cysll2 
(Figure ID and Additional file 1: SIl). The subsequent 
addition of pK displaced the anthraniloyl from Cysll2 
with HHQ formation as confirmed by mass spectrometry 
(see supporting information Text SI2). Strikingly, repeat- 
ing the experiments with inverted substrate order (pK 
first, then ACoA) resulted in less than half of HHQ forma- 
tion (Additional file 1: Figure SI2A) supportive for the 
preferential substrate order deducible from the kinetic 
analysis of the HHQ biosynthesis. However, the latter find- 
ing cannot be excluded to be at least in part due to sub- 
strate inhibition. The different plots of the kinetic data 
and the "inverted" SPR experiments sustain the idea that 
PqsD follows a ping-pong kinetic mechanism with ACoA 
as the first substrate (Figure IE). Based on these results a 
putative kinetic cycle for HHQ biosynthesis was set-up 
and different PqsD-ligand complexes chosen to simulate 
the single steps (Additional file 1: Figure. SI2B). 

Comparative analysis 

PqsD and KAS III enzymes are structurally similar, but 
they diverge in their biosynthesis pathways. Differences 
in the amino acid sequence of the active sites might 
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Figure 1 Ping-pong Icinetic mechianism of PqsD in HHQ biosyntiiesis. A) Lineweaver-Burk Plot and B) Hanes-Woolf Plot representing the 
activity of PqsD as a function of ACoA concentration in the presence of 60 (diamonds), 120 (black circles), 240 (triangles), 480 (open squares) and 
1000 |jM (open circles) BK. All data are mean values of four replicates. C) Kinetic parameters of PqsD. D) "Real-time" PqsD kinetics by SPR 
experiments: 1. ACoA injection; 2. buffer run - response units remain high due to anthranilate-CySii2; 3. (3K injection; 4. Buffer run - response line 
goes to zero due to HHQ formation and restored apo PqsD. E) Schematic view of the putative ping-pong kinetic mechanism of the HHQ 
biosynthesis (E - PqsD, Ecsj - PqsD with Cysl 12-ligated anthranilate). 



account for that. A BLAST search [34] of PqsD was 
performed against the Uniprot and the PDB database. 
The closest homologues were FabH2 of several patho- 
genic Burkholderia strains (-55% sequence identity) 
and FabHs of Gram-positive and Gram-negative bac- 
teria (-29-40% sequence identity). A multiple sequence 
alignment of PqsD with these enzymes was performed 
using Probcons [35] with Jalview [36] (Additional file 1: 
Figure. SI3-SI4). As seen in the sequence alignment all 
these enzymes, but in particular PqsD and Burkholderia 
FabH2, diverge in amino acid sequence in several re- 
gions (Figure 2 and Additional file 1: SI3). The domain 
that varies most in terms of sequence identity as well as 
in 3D-folding is the dimer-interface, comprising hairpin 
loop (hL; residues 185-220) and helixes H8-H9 (h8-9; 
residues 143-160). Further, differences are found in 
the ''substrate-loop' (sL) and the adenosine binding do- 
main (aBS). 

As evident in PqsD and FabHs 3D-structures the sLs 
border the active site and contribute in the formation of 



a central cavity in the dimer structures, which is 
surmounted by the hL helixes Hll and which can be 
filled by water molecules (Figure 2 and Additional file 1: 
SUB). Mutagenesis studies in FabH showed that resi- 
dues placed within this sL play an exceptional role for 
the substrate specificity [37]: bulky residues (Phe87; 
£. coli FabH - ecFabH) determine a clear preference for 
the short-chained acetyl-CoA, smaller residues (Thr87; 
Af. tuberculosis FabH - m^FabH) account for a long- 
chained substrate preference, such as lauroyl-CoA 
[37-40]. PqsD and Burkholderia FabH2 both produce 
HHQ and both share the same palindromic sequence 
(-S82PDHHDPS89-) in the sL. According to the se- 
quence alignment for PqsD and FabH2 we identified 
Asp87 as positional analogue of ecFabH Phe87 and 
m^FabH Thr87 (Figure 2); consequently, also, Asp87 
might be involved in the substrate recognition process. 

Mutagenesis studies on ecFabH showed that exchan- 
ging basic residues surrounding the primary funnel 
access to acid residues strongly affected CoA- and 
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Figure 2 Functionally important domains of PqsD. A) 3D structure of the PqsD dimer with functionally significant domains and the central 
cavity. PqsD is rendered as cartoon with flexible regions coloured in green. Water molecules present in the central cavity are shown as 
red spheres. B) Multiple sequence alignment of PqsD and homologous KAS-III enzymes. Abbreviations: Paer - Pseudomonas oeruginoso, 
Bamb - Burkholderia ambiforia, Ecoli - Escherichia coli, Aael - Aquifex aeolicus, Tther - Thermus thermophilus, Saur - Staphylococcus aureus, 
Efaec - Enterococcus faecalis, Sliv - Streptomyces lividans, Mtub - Mycobacterium tuberculosis. 



ACP-binding [41]. As shown in Figure 2 in PqsD the 
corresponding residues are negatively charged or neu- 
tral (i.e. Glu227, Glu269, Gln270; violet triangle), 
which is reflected on a tertiary structure level in a 
modified electrostatic potential surface (Additional file 1: 
Figure SIS). The access to the primary funnel of PqsD 
is surrounded by basic residues, in part not present 
in other FabHs, forming a large "cationic belt" (i.e. 
Arg36, Argl51, Argl53, Arg221, Arg223, Arg262). 
This basic network forms part of the CoA and ACP 
binding site [31]. 



Molecular dynamics simulations (A-F) of the main steps of 
the kinetic cycle 

MD simulations of nine PqsD complexes A-F were 
performed to elucidate the dynamic motions within the 
HHQ kinetic cycle and the binding modes of ACoA and 
pK. The volume variations of the primary funnels, sec- 
ondary channels and central cavity were tracked using 
the software fpocket2 [42] (Table 1). Finally, binding 
free energy differences (AGbind) for ACoA, (BK and 
HHQ were estimated following a single-trajectory ap- 
proach using MM-GBSA methodology (only enthalpy 
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Table 1 Internal pockets volume variations 



MD code Central cavity Primary funnel (chain A) Primary funnel (chain B) Secondary channel (chain A) Secondary channel (chain B) 
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Summary of the volume profiles (in A^) over the time for central cavity, primary funnels and secondary channels of chain A and B for each simulation B-F tracked 
with fpocket2. The trend progression of the volume is summarized as follows: decrease (<), modest decrease (<), constant (=), modest increase (>), increase (>). 
The time-dependent volume-profiles are added in supplementary information as Figure. SI10. 

MD codes: B - PqsD apoform, CI - PqsD-ACoA in chain B, C2 - PqsD with ACoA in the primary funnel of chain A and B, D - PqsD with CysTi2-bound anthranilate 
(PqsD-CSJ), El - PqsD-CSJ with pK in the primary funnel of chain B, E2a - PqsD-CSJ with pK in the secondary channel of chain A, E2b - PqsD-CSJ with pK in the 
secondary funnel of chain A and B, F - PqsD with HHQ in the primary funnel of chain B. 



was computed; Table 2) and performed an MM-GBSA 
binding free energy decomposition analysis to high- 
light relevant residues for ligand binding. As two bind- 
ing sites exist in the PqsD dimer we also inquired, 
whether synergistic effects resulted for AGbind compar- 
ing the simulations with two ACoAs (C2) and |3Ks 



Table 2 Estimated binding free energies (AGbind) using 
MM-GB/PBSA methods 



MD 
code 




GBSA 


E_GAS 


EGB AGbind STD 


CI 




PqsD + ACoA (chain B) 




-308 


196 -111 14 


C2 




PqsD + ACoA (chain A) 




-32 


-50 -82 8 






PqsD + ACoA (chain B) 




25 


-88 -63 7 


El 


PqsD-CSJ + single (JK, primary funnel (chain A) 




-44 


17 -27 2.5 
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PqsD-CSJ + single (JK, secondary channel (chain A) 




15 


-42 -27 6 
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PqsD-CSJ + PK in secondary channel (chain A) 




-7 


-37 -44 5 




PqsD-CSJ + PK in secondary channel (chain B) 
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-32 -26 6 


F 




PqsD + HHQ (chain B) 




-54 


14 -40 3 



STD. - standard deviation; E_GAS - binding energy in vacuum, EGB - binding 
energy in implicit solvent computed with GB method. 
AGbind (kcal/mol) are computed for ACoA, pK and HHQ in the different PqsD 
complexes describing the kinetic cycle: CI - PqsD with ACoA in chain B, C2 - 
PqsD with ACoA in the primary funnel of chain A and B, D - PqsD with Cysn2- 
bound anthranilate (PqsD-CSJ), El - PqsD-CSJ with pK in the primary funnel of 
chain B, E2a - PqsD-CSJ with pK in the secondary channel of chain A, E2b - 
PqsD-CSJ with pK in the secondary funnel of chain A and B, F - PqsD with 
HHQ in the primary funnel of chain B. 



(E2b) to their single ligand counterparts (CI and E2a 
respectively). 

Flexibility of PqsD as monomer or dimer. We investi- 
gated the flexibility of PqsD submitting monomer A of the 
apoform structure 3H76 to the HingeProt [43] and the 
HingeMaster [44] web-servers. Both predictions indicate 
the existence of hinges at the N- (G185/T195) and the 
C-terminal end (M220/G222) of the hL capable of large 
conformational changes (Additional file 1: Figure SI6A-C; 
morphing the transition from closed to open hL con- 
former using the Yale Morph Server [45]; Additional file 2: 
Movie SI: video representing hinge opening and closing). 
A more realistic measure of the PqsD flexibility is obtained 
from the monomer (A) and, especially, the dimer PqsD 
apoform MD simulations (B): in the monomer MD 
very large fluctuations are observed for the hL (Figure 3; 
green line), which collapses closing the secondary channel 
(Additional file 1: Figure SI7A). In the dimer MD, on con- 
trary, hL and sL of both chains sustain each other with 
only the C-terminal |3-sheet of the hL moving out thus ex- 
posing the secondary channel access. The hinge predic- 
tions and the motions seen for these MD simulations 
evidence the flexible nature of the hL and suggest an im- 
portant role for it in PqsD functionality. 

Accounting for all the MD simulations (B-F) a rather 
broad conformational ensemble was gathered. Thereby, 
most of the fluctuations were located in the upper third 
of PqsD (Additional file 1: Figure SI7B) in regions also 
found disordered in several FabH crystal structures. This 
implies that conformational rearrangements in any of 
them can significantly impact on the functioning of the 
enzymes. The most flexible domain in all the MD simu- 
lations is the dimer interface including hL and h8-9 
(Additional file 1: Figure SIS). 

The refolding of the two hL helixes Hll and of h8-9 
affected the dimer interface morphology (yellow dotted 
circle - cavities at 0 ns, yellow full circle - cavities at 
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Figure 3 Residue-dependent RMS fluctuations of the apoform MD simulations with PqsD in its monomeric and dimer state. The 

fluctuations of the residues of PqsD in its monomer form (MD code A) are coloured in green, whereas for the PqsD in its dimer state (MD code 
B) chain A and B are shown in blue and red respectively. Domains putatively involved in the catalysis as identified by comparative analysis with 
homologous KAS-III enzymes are labelled with their abbreviations (dotted orange box indicates helix HI 4): adenosine binding site - aBS, 
palindromic "substrate-loop" - sL, helix H8-H9 loop - h8-9, hairpin-loop - hL, helix HI 2 - hi 2, "oxyanion-loop" - oL. Notably, also for all other MD 
simulations the largest RMS fluctuations were found in these areas (see also Additional file 1: Figure SIS). 



30 ns; Figure 4) and the geometry of the arginines sur- 
rounding the primary funnel access ("cationic belt"; 
Additional file 1: Figure SI9). Thereby, deep channels are 
formed towards the central cavity (for apoform B, 
single-ACoA CI and the CSJ-PqsD D simulations) or 
opening out into the secondary channel (for simulations 
with pK in the secondary channel E2a-E2b; Figure 4). 

In the apoform trajectory B the enzyme floats between 
different conformers. On contrary, in presence of the 
different substrates the equilibrium is selectively shifted 
toward one preferred structure. This is reflected in the 




Figure 4 Conformational changes at dimer Interface In the 
dual-^K MD simulation E2b. Top view of the dimer interface at 
30 ns with the secondary channel solvent-exposed (dotted yellow 
eclipse) and the (3K molecules shown as green surfaces. The dashed 
green line indicates a cleft which exists between the two HI 1 
helixes at 0 ns, but that disappears over the simulation run. The full 
yellow circles indicate additional grooves formed at the interface. 
The cyan arrows indicate the primary funnel accesses for 
both chains. 



volume variations of internal cavities and channels 
(Table 1 and Additional file 1: Figure SIIO). In particular, 
tracking the central cavity volume of the MD simula- 
tions B-F resulted in a sinusoidal, "heart-beat"-like 
volume-profile (Table 1): in the apoform simulation (B) 
the volume is reduced; in the single-ACoA complex 
(CI) the cavity volume increases; with two ACoA mole- 
cules (C2) as well as for the CSJ-PqsD (D) complex the 
cavity volume decreases. Finally, in the single- |3K E2a 
and, in particular, in the dual-pK simulations E2b the 
volume increased. The formation of channels from the 
central cavity into the oxyanion site (close to Asp87 and 
Ser317), sideways out to the surface, or up to the dimer 
interface suggests that water molecules might move 
along these channels depending on the catalytic needs. 

The fact that ligand-induced structural changes are 
observed in the very same regions in presence of diverse 
ligands makes us confident that a sufficiently large con- 
formational ensemble has been gathered by the MD 
simulations to represent a good starting point for 
structure-based drug design and virtual screening. A 
good example is represented by the motions of Phe218, 
which is located on the C-terminal |3-sheet of the hL of 
each monomer (Additional file 1: Figure Sill): in the 
MDs with ACoA CI and C2, with CSJ D, with |3K in the 
primary funnel El, and with HHQ F (chain B) it stays 
turned towards the catalytic centre occluding the access 
from the catalytic site to secondary channel and dimer 
interface. On contrary, when pK is in the secondary 
channel (E2b and chain A of E2a) Phe218 rotated out- 
wards increasing the distance between the centre of 
mass of Phe218 and the Ca of Cysll2 (Additional file 1: 
Figure Sill); the secondary channel opened out to the 
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dimer interface and the catalytic centre was enlarged. 
When no ligand is present in the catalytic centre (e.g. in 
the apoform MD B and in chain A of F), an intermediate 
position can be observed indicative for the conform- 
ational flexibility of this area. 

ACoA progression into the catalytic site 

The simulations with ACoA supplemented well the 
PqsD crystallographic data. In particular for the single- 
ACoA simulations CI ACoA moved deep into the pri- 
mary funnel ending in a pose where the anthranilic ring 
overlaps well with the Cysii2-bound anthranilate of the 
PDB-ID 3H77 structure (Figure 5 and SI12). Strikingly, 
only Cysll2, His257, and Asn287 of the monomer 
hosting ACoA (chain B) assumed an ideal geometry for 
catalysis (e.g. Cysll2-His257 ~4 A; initial/final distances 
in Additional file 1: Figure SI12B). Arg36, Argl53, and 




Figure 5 Binding mode of ACoA. A) Final poses of the 
single-ACoA simulation CI (cyan; 37 ns) and of the dual-ACoA 
simulation C2 (green; 31 ns) are superimposed with the PDB 
structure 3H77 (ACoA orange - starting position; anthranilate-CySii2 
in magenta - overlaps with ACoAci). Phe32, Argl53, Asnl54, 
Phe218, His257, Asn287 and Ser317 of the single-ACoA simulation at 
37 ns are shown as dark cyan sticks. B) Decomposed energy 
contributions per residue (at least for one MD >0.5 kcal/mol) 
determined by MM-GBSA methods for the MD simulations CI (cyan) 
and C2 (chain A - red; chain B - black). 



Arg223 showed the highest binding energy contributions 
due to interactions with the phosphate groups of the 
pantetheine-linker. As obvious from the decomposed en- 
ergy contributions (Figure 5B) also Thr28, Phe32, Cysll2, 
Asnl54, Leul55, Ilel57, Leul58, Leul93, Phe218, Met220, 
Gly222, Met225, Phe226, His257, Pro259, Asn260, Ile263, 
Asn287, Ala289, Tyr315, and Ser317 contribute to ACoA 
binding. The optimal position reached by ACoA in the 
simulation CI is also reflected in a more favourable esti- 
mated AGbind (-111 kcal/mol; chain B) than for the two 
ACoA molecules of the dual-ACoA simulation C2 
(-82 kcal/mol, chain A; -63 kcal/mol, chain B). Strikingly, 
plotting AGbind-versus-time shows a progressively more 
favourable AGbind for CI (Additional file 1: Figure SI13), 
whereas for the two ACoAs in C2 the AGbind remain ra- 
ther constant indicative for an equilibrated system, in 
which, however, the ACoAs in C2 do not reach a 
catalysis-like position (see distances in Additional file 1: 
Figure SI12C). This hints at a sequential ACoA entrance/ 
release mechanism, where the empty monomer refolds 
while binding of the first ACoA. 

Where and how does |3K bind in PqsD? 

Our SPR experiments showed that using free pK acid as 
second substrate yields higher HHQ formation than when 
it is added first. The active site of PqsD PDB structures 
seems inadequate to handle pK-(ACP) binding and pK- 
incorporation [23] . However, in PqsD a secondary channel 
exists similar to that of mtFabH, which can host long- 
chained p-keto-acids (i.e. lauroyl-CoA) [40]. The access to 
this channel is lined by polar residues (Argl45, Thrl95, 
Ser317, and Asp87 of the second monomer), whereas its 
bottom part is rather hydrophobic (Leu81, Leul42, 
Leul55, Leul59, Leul93, Metl94, Phe218 and Met220). 
Further, this secondary channel borders the central cavity 
and the dimer interface with the ion-pair Asp87-Argl45 
(of the other chain) and Phe218 acting as barrier respect- 
ively. As PqsD is clearly capable of HHQ biosynthesis con- 
formational changes are expected, in analogy to mtFabH 
[46], that allow HHQ formation and release. 

We followed two approaches to determine the most 
plausible access path and binding mode of pK. In the first, 
we docked pK into the primary funnel of anthranilate- 
ligated PqsD. Two main orientations resulted: 1) with the 
polar head of pK pointing to the catalytic triad (Additional 
file 1: Figure SI14A); and 2) turned 180° with the carbox- 
ylic group interacting with the Arg of the cationic belt. 
The second binding mode was not supportive for any of 
the catalytic mechanisms [10,17]. 

The second approach was based on the substrate size 
similarity and the structural homology between PqsD 
and m^FabH. We superimposed the Cysll2-anthranilate 
PqsD (PDB-ID 3H77) and mtFabH co-crystallized with 
dodecyl-CoA (PDB-ID 1U6S) [46] and then transposed 
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and modelled the 3-oxo-undecanoyl from the mtFabH 
structure into pK in the secondary channel of PqsD 
(Additional file 1: Figure SI14B). Active site refinement 
with the LigX module of MOE [47] with restrained |3K 
resulted in an energy optimized complex with small Ca 
RMSD (~ 1 A) compared to the starting complex. This 
second binding mode of |3K looks similar to that 
obtained by Bera et al [23] for decanethiol by superim- 
posing PqsD with the m/:FabH structure PDB-ID 2Q01 
where the decanethiol is covalently attached to the cata- 
lytic Cys. 

The MD simulations El-E2b with pK either in the pri- 
mary or in the secondary channels were supportive to 
understand the putative ping-pong kinetic cycle. 

MD simulation CSJ-PqsD with p/C in the primary fun- 
nel of chain B (El; 34 ns): Only small conformational 
changes are observed in this MD simulation for pK, 
which points the 3-oxo-p-keto moiety towards the 
catalytic triad (Additional file 1: Figure SI15). Thereby, 
the carboxylate head is trapped in a hydrogen-bond 
network with His257 and Asn287 holding the p carbon 
of pK close to the CSJ-sulphur as shown by their 
favourable energy contributions (Additional file 1: 
Figure SI15B). Additionally, van der Waals interactions 
are formed between pK and Leul93, Met220, Met225, 
Phe226, Pro259, Ile263, and Tyr315 (Additional file 1: 
Figure SI15). 



MD simulation CSJ-PqsD with pK in the secondary chan- 
nel of chain A (E2a; 30 ns): After 30 ns the carboxylic moi- 
ety of pK protruded into the catalytic centre within 4.5 A 
of the amine of CSJ, the hydroxy group of Ser317 and of 
the oxyanion site formed by Glnlll and Ser317. Further, 
pK forms a hydrogen bond with Thrl95 as shown in 
Figure 6. Phe218 was pushed outwards by pK in chain A 
(Figure 6A and Additional file 1: Figure Sill), while in the 
empty chain B Phe218 points towards the catalytic centre 
(Additional file 1: Figure Sill), reducing the primary fun- 
nel volume of chain B (Additional file 1: Figure SI 10). 

MD simulation CSJ-PqsD with p/C in both secondary 
channels (E2b; 30 ns): As seen comparing the single- and 
the dual-ACoA MD simulations CI and C2 (Additional 
file 1: Figure SI 17 and SI 18) simulating two pK molecules 
simultaneously resulted in a more stable trajectory 
(smaller RMSD as for E2a; Additional file 1: Figure SI21 
and SI22). Still, the modelled presence of pKs in the sec- 
ondary channels induce structural changes in the 
dimer-interface and in neighbouring domains (hL, h8- 
9, hl2 and oL of one monomer and sL of the other 
monomer; Additional file 1: Figure SI22). Two long 
U-shaped tunnels are formed, which connect primary 
funnel access, catalytic site (with CSJ and pK in a pro- 
ductive pose) and secondary channel access at the dimer 
interface (Figure 7). The two pK are in a similar or even 
more optimal pose (see pK in chain A) as described for 




Figure 6 Binding mode of PK in the secondary channel (MD simulations E2a and E2b). A) Superimposed snapshots of (3K binding to PqsD 
at 0 and 30 ns. (3K at 0 and 30 ns is sliown as purple sticl<s and as green ball and sticks respectively. Important residues are shown as cyan lines 
at 30 ns and as orange lines at 0 ns. The red arrows indicate the conformational shift from the initial (yellow) to the final (blue) conformation of 
PqsD in the MD simulation E2b; Phe218 flips out in presence of PK. B) Schematic representation of PK in the secondary channel. Polar amino 
acids are illustrated in purple and hydrophobic amino acids in green circles. Hydrogen bonds and CH-pi interactions are shown as green arrows 
and dotted lines. C) Decomposed energy contributions per residue (at least for one MD >0.5 kcal/mol) determined by MM-GBSA methods for the 
MD simulations E2a (cyan) and E2b (chain A - red; chain B - black). 
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the single- |3K MD E2a. Phe218 is turned outside 
(Figure 6A and Sill) exposing the secondary channels 
of chain A and B to the dimer-interface (Figure 4). 
Hydrogen bonds are formed between the carboxylic 
moiety of pK and the amine of CSJ and the hydroxy 
group of Ser317. Further, a hydrogen bond is formed 
between the p-keto carbonyl and Thrl95 holding pK 
close to the "substrate-loop" (sL) and in particular to 
Pro88 and Asp87 (see binding energy contributions; 
Figure 6C). Finally, Van der Waals interactions are 
established between pK and Leu81, Pro88, Leul55, 
Leul59, Leul93, Metl94 and Phe218. 

Catalytic mechanism and |3K binding modes 

In this study the free acid form of pK was used to test 
HHQ biosynthesis. Two aspects have thus to be consid- 
ered: firstly, the free pK-acid is unlikely to exist in bacteria 
and its /<d value in the millimolar range determined by 
SPR and enzyme kinetics does not fit with the substantial 
HHQ production in aeruginosa. Secondly, inverting the 
substrate order in either the SPR or the enzyme kinetic 
studies still yields HHQ formation, although substantially 
decreased (see Additional file 1: Figure SI2A). 

The final pose identified in the simulation El with pK 
in the primary funnel hints at a potential pK-activation 
via decarboxylation, which could occur either spontan- 
eously or facilitated by His257 and Asn287 in analogy 
to other KAS-III enzymes.24 This pose of pK fits well 
with the decarboxylative Claisen reaction (Figure 8E). 



However, both of the above-mentioned aspects argue 
against this binding mode as the "m vivo' one: the 
orientation of pK cannot match that of any thioester- 
bound form needed for its delivery (Figure 8A), and, in 
addition, no access to the catalytic site for ACoA as sec- 
ond substrate is possible as long as pK is placed in the 
primary channel as seen for El (Figure 8B). 

Given the unlikeliness of free pK acid as "m vivo' sub- 
strate other events must occur in HHQ biosynthesis, 
such as involvement of PqsB and PqsC [31] or ACP- 
thioester binding, which could reduce the activation en- 
ergy or favour the kinetics. Holding true the ACP- 
delivered pK, we notice that the binding site topology 
and final orientation of pK in the secondary channel 
(MDs E2a-E2b) does not disturb the access of ACoA 
(Figure 8C) nor it precludes thioester-delivered pK 
(Figure 8D). The interactions with CSJ112, Thrl95, and 
Ser317 found for E2a-E2b (see energy contributions in 
Figure 6C) fit well with the imine/enamine mechanism 
proposed by Diggle et al. [10] (Figure 8F). An equal or 
even more favourable AGbind is estimated for the two pK 
in the dual-pK simulation E2b (-44 kcal/mol - chain A, - 
26 kcal/mol - chain B) as compared to AGbind of the 
single-pK (-27 kcal/mol in El; -27 kcal/mol in E2a) 
(Table 2). This hints at a positive synergistic effect of the 
simultaneous presence of two pK molecules in the second- 
ary channels, which might facilitate conformational 
changes necessary to accommodate pK in an energetically 
favoured pose. In this regard, it is helpful to remember 




Figure 7 U-shaped channel network formed in presence of PK in the secondary channel. The U-shaped channel connects primary funnel 
access with the catalytic centre (large unfilled arrow), the secondary channel (thin arrow) and the dimer interface (large dotted arrows). In 
presence of (3K a tunnel is formed between Cysii2-bound anthranilate (CSJ) and the central cavity (dashed thin arrow) potentially used for water 
supply (needed for catalysis) from the central cavity to the oxyanion site close to P88'. pK is shown as green ball and sticks, whereas residues 
shown to interact with pK (CSJ, T195, F218, RMS, D87', P88') as cyan sticks. For clarity the domains are labelled with their abbreviations. The 
solvent accessible surface is color-coded as follows: blue - positive electrostatic potential (+25 kcal/mol); red - negative electrostatic potential 
(-2S kcal/mol), white - neutral. 
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Figure 8 Binding modes of PK and putative catalytic mechanisms of PK-incorporation. A) (3K is placed in tlie primary funnel (blue tube) as 
in the MD simulation El: no space is left to host the thioester carrier CoA/ACP. B) The access of ACoA into the primary funnel is hindered with 
|3K in the pose as in A. C) (3K is in the secondary channel (violet box) as in the MD simulations E2a and E2b: in the primary funnel space is 
available to accommodate the thioester carrier CoA/ACP. D) ACoA can access the primary funnel with (3K in this pose and transfer anthranilate to 
Cysl 12. E) The Claisen condensation fits best as mechanism for PK incorporation for (3K in the primary funnel. F) The imine/enamine formation is 
the most plausible mechanism for (3K incorporation when (3K is in the secondary channel. 



that in the dual-pK simulation a dynamic channel-system 
is formed, which enlarges the catalytic centre (Additional 
file 1: Figure SIIO), thus accounting for a putative intramo- 
lecular cyclization in HHQ formation. Although no con- 
clusive data exists all the above aspects let us favour the 
secondary channel binding mode as the more probable. 



due to the rooted ability of S317F to form hydrogen 
bonds, as suggested by the consistent binding energy 
contributions of Ser317 in most of the MD simulations 
with both ACoA and pK (Figure 5 and 6), or because of 
the primary funnel being occluded by the newly intro- 
duced aromatic ring. 



Ser317Phe PqsD mutant 

In the MD simulations Ser317 is involved in hydrogen 
bonds with the carbonyl-group of ACoA (CI and C2) 
and with the carbonyl- or the carboxylic-moiety of |3K 
(E2a and E2b). To verify its importance we exchanged 
Ser317 by site-directed mutagenesis into Phe and deter- 
mined the catalytic activity by detection of the formed 
HHQ using UHPLC-MS/MS. Under assay conditions 
(0.1 (iM enzyme, 5 (iM ACoA, 70 (iM |3K) the mutant 
produced less factor 700 HHQ compared to the wild 
type. Tenfold higher enzyme concentration of the S317F 
did not result in an increase in HHQ production indicat- 
ing a complete abolishment of the catalytic functionality. 

In order to understand at which stage of the kinetics 
the S317F mutant disrupts the biosynthesis we per- 
formed SPR experiments using S317F PqsD. Strikingly 
no anthranilate transfer to Cysl 12 was detected (Fig- 
ure 9) suggesting that substitution of Ser317 with Phe 
already disrupts the first step of the kinetics. The close 
vicinity of Ser317 to the catalytic site, however, makes it 
questionable whether the very weak enzyme activity is 



Conclusions 

In this work we elucidated some of the molecular bases 
of HHQ biosynthesis. A putative ping-pong kinetic 
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Figure 9 SPR sensograms of ACoA. When ACoA is injected with 


PqsD wild-type fixed on the chip (black) an increase in the response 


units is observed corresponding to the anthranilate transfer to 


Cysl 12. Differently, when S317F PqsD mutant (red) is present on the 


chip, no increase is observed. 
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mechanism was determined by enzyme kinetics experi- 
ments, which was substantiated by the preferred sub- 
strate order (ACoA prior to |3K) assessed by SPR. These 
data were useful to set-up and analyse the MD simula- 
tions, which aimed to unmask dynamic motions 
governing PqsD functionality. MD simulations reveal a 
kinetic-step dependent adaptation of PqsD to the di- 
verse ligands, a favoured binding mode for |3K in the 
secondary channel, as well as an arginine-assisted pro- 
gression of ACoA towards the catalytic site. Also, 
Ser317 was identified as an important binder for ACoA 
and |3K, at least in part confirmed by the inactive S317F 
PqsD mutant. 

In general, SPR studies with PqsD mutants will be 
helpful to determine the binding mode of inhibitors. Fi- 
nally, the conformational ensemble retained from the 
MD simulations will serve as valuable starting point for 
structure-based design of PqsD inhibitors. 

Methods 

Production and purification of recombinant PqsD in E. coli 

The overexpression and purification of PqsD as an N- 
terminal Hise-tagged fusion protein in E, coli BL21 (XDE3) 
using the vector pET28b(+) -pqsD was performed as de- 
scribed by Pistorius et al. [17]. To remove the Hise-tag the 
protein was subjected to thrombin cleavage performed at 
16°C for 16 h in a 50 mM Tris-HCl buffer, pH 8.0, 
containing 150 mM NaCl, 1 mM 2-ME, 2.5 mM CaC^ 
and 1 unit thrombin per mg protein followed by a second 
passage through the His Trap HP 5 column. The protein 
was frozen in liquid nitrogen. 

Preparation of S317F PqsD mutant 

S317F PqsD mutant was generated using the QuikChange 
Site-Directed Mutagenesis Kit (Stratagene, La Jolla, CA) 
according to the manufacturers instructions using the 
pET28b(+)/j9<5^5£) plasmid as template and the mutagenesis 
primers 5' GCTGGTCCTGACCTACGGTTTTGGCGC 
GACCTGGGGCG 3' and 5'CGCCCCAGGTCGCGCC 
AAAACCGTAGGTCAGGACCAGC 3'. Plasmid DNA 
was purified using the GenEluteTM HP Plasmid Miniprep 
Kit ( Sigma- Aldrich, St. Louis, MO) and sequenced to con- 
firm the site-directed mutation. 

Enzyme kinetic analysis 

The PqsD catalysed formation of HHQ was analysed 
using a UHPLC-MS/MS based assay performed in 96- 
well microtiter plates (Greiner) using the method of 
Pistorius et al. [17] with some modifications. First, the 
purified enzyme PqsD (0.8 (iM; in 50 mM MOPS, 
pH 7.0, 0.016% (v/v) Triton X-100) was preincubated 
without substrates for 5 min at 37°C. The reaction buffer 
(15 (iL; 50 mM MOPS, pH 7.0) and the substrates ACoA 
(20 (iL; concentrations: 2-40 (iM) and EK (20 (iL; 



concentrations: 240-4000 (iM) were added. The reaction 
was started by the addition of preincubated enzyme 
(25 (iL; 0.8 (iM) resulting in a total reaction volume of 
80 (iL with the following final concentrations: PqsD 
0.25 (iM, ACoA 0.5 - 10.0 (iM and £K 60-1000 (iM, Tri- 
ton X-100 0.005%, methanol 2%. The reaction was 
stopped after 4 min at 37°C by adding 80 (iL of metha- 
nol containing 1 (iM of the internal standard amitriptyl- 
ine. For each sample, the reactions were performed in 
quadruplicate. HHQ-formation was detected using 
UHPLC-MS/MS according to method in supplementary 
information (Additional file 1: Text SIl). Data were plot- 
ted and analysed using GraphPad Prism 5 software. 

Synthesis of Anthraniloy IS -Coenzyme A thioester 
{ACoA), ACoA was synthesized from isatoic anhydride 
and coenzyme A (CoA) as described by Simon and 
Shemin [48]. 

Synthesis of ethyl 3'Oxodecanoate, S-oxodecanoic acid 
(|3 ketodecanoic acid), and of HHQ {2'heptylquinolin'4 
{lH)'One), Synthesis as described by Lu et al. [49]. 

Surface Plasmon resonance (SPR) 

SPR binding studies were performed using a Reichert 
SR7000DC optical biosensors instrument (Reichert 
Technologie, Depew, NY 14043 USA). HClOOOm sen- 
sor chips were purchased from Xantec Analytics 
(Diisseldorf, Germany). 

Immobilization of HePqsD or HePqsD S317F 

Overexpression and purification of PqsD or of the S317F 
PqsD mutant was performed as previously described 
[32]. HePqsD or H^PqsD S317F was immobilized on 
HClOOOm sensor chip using standard amine coupling 
chemistry at 25°C analogous to the manufacturers in- 
struction. H^PqsD or H^PqsD S317F was diluted into 
10 mM sodium acetate (pH 4.5) to a concentration of 
100 [ig/mL and coupled to the surfaces with densities 
between 15000 and 20000 RU. 

Catalytic activity of PqsD 

ACoA was diluted into running buffer to 10 (iM. 
|3-ketodecanoate (|3K) was dissolved in methanol to a 
10 mM stock solution and diluted into running buffer to 
20 [iM, In the first experiment ACoA was injected for 
5 min association and 10 min dissociation time, followed 
by a 20 min injection of pK. In the second experiment 
|3K was added to the running buffer (100 [iM in 50 mM 
Tris, pH 8.0, 150 mM NaCl, 0.1% Triton X-100). The 
ACoA-injection (10 \iM for 5 min) followed when the 
binding site was saturated, indicated by a stable baseline. 
In all experiments the flow through was collected. 
After ACoA addition CoA emission was detected 
using UHPLC-MS/MS as reported elsewhere [32]. 
After addition of the second substrate (experiment 
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dependent) the flow through was extracted with 1 mL 
CHCI3, evaporated and diluted in 50 (iL MeOH. HHQ 
formation was detected using UHPLC-MS/MS (see 
Additional file 1: Text SIl). 

SPR study with S317F PqsD mutant 

For comparison of ACoA binding to PqsD wild- type 
and to S317F PqsD mutant, the study was performed as 
in the first experiment described above. 

Computational methods 

Docking of p/C and HHQ, |3K and HHQ were docked 
each 50 times with GOLDv5.0 [50] using the GOLD 
SCORE function [51]. The docking site was defined in- 
cluding all residues within 9 A of ACoA found in PDB 
structure 3H77. Covalent ligated PqsD and apoform 
PqsD were used to dock |3K and HHQ, respectively. The 
default GOLD parameters were used. 

Simulation protocol 

System setup was performed as follows for all simulated 
systems. Atomic coordinates were taken from PDB-ID 
3H76 for MD simulations A, B, CI, C2 and F, and from 
3H77 for D, El, E2a and E2b. Water molecules and ions 
present in the crystals were removed. The protonation 
states were determined at pH 7.4 with the Protonate3D 
module of MOE and the enzyme complexes then mini- 
mized using the MOE module ligX [47]. The solvated 
systems were set up using the AMBERll [52] program 
xLeap with AMBER99SB force field [53]. A 9 A pad of 
TIP3P waters was added to solvate each system as octa- 
hedral box. Neutralizing counter ions were added to 
each system. 

Parameters for ACoA, pK, HHQ and the Cysll2- 
bound anthranilate (CSJ) were determined using the 
sqm routine of AmberToolsl.5 [52]. For each ligand 
AMl-BCC charges were computed. For ACoA and pK a 
net charge of -3 and -1 was set, respectively. For CSJ, 
which was taken from the PDB-ID 3H77, the Cys pa- 
rameters were used as starting point followed by 
parameterization with the sqm routine. 

MD simulations were performed with the parallelized 
PMEMD module of AMBER 11. The starting PqsD com- 
plexes were minimized and equilibrated with the back- 
bone atoms restricted by harmonic restraints of initially 
10 kcal mol AT^ and then progressively reduced to 5, 
2, 1 and 0 kcal/mol. The systems were heated to 300 K 
in the canonical NVT ensemble (constant number of 
particles, N; constant volume, V; constant temperature, 
T) using a Langevin thermostat, with a collision fre- 
quency of 3.0 ps"^ A"^. Production runs were then made 
for 30-37 ns duration in the NPT ensemble at 300 K. 
As with the heating, the temperature was controlled 
with a Langevin thermostat (but with a 1.0 ps"^ collision 



frequency). The time step used for all stages was 2 fs 
and all hydrogen atoms were constrained using the 
SHAKE algorithm [54]. Long-range electrostatics were 
included on every step using the Particle Mesh Ewald 
algorithm with a 4th order B-spline interpolation [55]. 

MD analysis 

B-factor, distances and RMSD time series were calcu- 
lated using the cpptraj analysis tool of the AmberTools 
1.5 package. 

Structures were sampled at 20 ps intervals. RMSD 
values were calculated for the enzyme, but also for each 
of the following regions individually, facilitating the in- 
terpretation of the data: adenosine binding site (aBS), 
"substrate-loop" (sL), helix8-9 (h8-9), hairpin loop (hL), 
helix 12 (hi 2), oxyanion loop (oL). 

Volume variations of internal cavities 

The volume variations (A) over the time (ns) for five 
cavities in PqsD were tracked for all PqsD dimer MD 
simulations using the software package fpocket2 [42]. 
All MD snapshots were superimposed using the C^ of 
the apoform PDB-ID 3H76 as reference structure. De- 
fault parameters for the identification of small cavities 
and channels were used. All plots are added in supple- 
mentary materials. 

MM-GBSA calculations 

Binding energies AGbind for ACoA, pK and HHQ were 
estimated by conventional MM-GBSA methods (Mo- 
lecular Mechanic - General Born Surface Area) [56] 
using snapshots of the simulations sampled every 30 ps. 
Energy decomposition analyses with Generalized Born 
solvent were performed with per-residue decomposition 
and 1-4 interactions added to the electrostatic and Van 
der Waals terms (idecomp = 2). More details to the 
MM-GBSA method are in Additional file 1: Text_SI2. 

A detailed description of the MD simulations A-F, cavity 
volume variation- ver5M5- time profiles (Additional file 1: 
Figure SI7), time-dependent AGbind profiles (Additional 
file 1: Figure SIIO), and RMSD plots (Additional file 1: 
Figure SI13-SI20) are added as additional data. The RMSD 
analyses of the different regions helped to visualize mo- 
tions in the MD simulations not visible from the all-atom 
RMSD plots for the entire protein and complemented the 
amino acid fluctuation analyses shown in Figure 3. 

Figures and plots 

Plots were made with Origin 9 or Excel, while figures 
with PyMOL [57] or MOE [47]. 
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Additional files 



Additional file 1: Supplemental methods, figures and references. The 

file is in PDF format. It includes: Text SI1 - UHPLC-MS/MS analysis of HHQ. 
Text SI2 - MM-GBSA theory. Text SI3 - Description of molecular dynamics 
simulations A-F. Supplementary Information Figures. Figure SI1 - Structural 
peculiarities of PqsD. Figure SI2 - HHQ biosynthesis follows a ping-pong 
mechanism. Figure SI3 - Multiple sequence alignment of PqsD and 
KAS III enzymes. Figure SI4 - Average distance tree from the ProconsWS 
alignment. Figure SIS - Electrostatic potential of PqsD (A) and £ coli FabH (B). 
Figure SI6 - Predicted hinge regions in PqsD. Figure SI7 - PqsD flexibility in 
the single monomer MD simulations A. Figure SIS - Residue-dependent RMS 
fluctuations for the MD simulations A-F. Figure SI9 - Conformational changes 
of the cationic belt. Figure SI10 - Time-dependent volume variations of 
internal cavities. Figure Sill - Time-dependent distance variation between 
Phe218 and Cysl 12. Figure SI12 - Progression of ACoA in the single-ACoA 
MD simulation CI. Figure SI13 - Time-dependent variation of the 
estimated binding free energy. Figure SIM - Where does (3K bind in 
PqsD? Figure SI15 - Binding mode of pK in the MD simulation El. 
Figure SI16-S23 - Trajectory analysis of the MD simulations B-F. 
Supplementary information References. 

Additional file 2: Movie 51. The morphing from the closed to the 
open hairpin-loop (hL) conformation is showed as result of the 
YaleMorphServer. The file is in avi format. 
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